Numerical and Experimental Analysis of Stress–Strain Characteristics in DP 600 and TRIP 400/700 Steel Sheets

The body constitutes the largest proportion of the total vehicle weight. Recently, increasing efforts have been made towards reducing its weight and improving its crashworthiness. By reducing its weight, fuel consumption will be reduced, and this will also translate into lower CO2 emissions. In terms of safety, vehicle body components use high strength steel which can absorb a substantial amount of impact energy. The present study pays attention to DP 600 and TRIP 400/700 stress–strain characteristics at quasi-static strain rates. The stress–strain characteristics of absorption capacity, stiffness, and deformation resistance force were investigated experimentally by tensile tests, three-point bending tests, and numerical simulations. The results indicate the potential for increasing the absorption capacity, stiffness, and deformation resistance force of the vehicle body’s deformable steel components. The present study verified the possibility of replacing physical testing with numerical simulation. A reasonably satisfactory agreement between the experimentally determined stress–strain characteristics and the numerical simulation was achieved, which can reduce the development time of deformable vehicle body components, reduce costs and optimize the selection of materials. The results extend the state of knowledge on the deformation characteristics of high-strength materials and contribute to the optimization of body components in terms of passive safety and weight.


Introduction
In the automotive industry, major efforts are made to develop coherent and balanced car concepts combining safety, emissions and affordability requirements.Such a combination of factors is required to be achievable through realistic, technologically and affordably viable production methods while avoiding an excessive burden on the environment.One of the factors that largely determines the success or failure of a given car concept is its safety.In terms of passive safety, a crucial role belongs to the car body, which consists of a front crumple zone, a safe zone for passengers (safety cage) and a rear crumple zone.In the event of a collision, the deformation structure of the vehicle body's front and rear sections is designed to absorb as much kinetic energy as possible so that the permissible biomechanical limit of overloading the human body is not exceeded in the event of a collision.The deformation of the structure at the front and rear of the vehicle is controlled.The forces generated by the impact in the bumper cross member are distributed to the longitudinal beams.From the longitudinal beams, the stress is transferred to the tunnel, to the sills, to the floor, to the door structure, and through the door pillar to the roof structure beams.In doing so, the longitudinal beams are subjected to a compressive force and the transverse beams and stiffeners are subjected to a bending force.The consequences of a side impact accident are often more severe than those of a frontal impact because the deformation structure of this part of the vehicle has to absorb the impact energy over a shorter path than in a frontal impact [1][2][3][4][5].
The absorption capacity (EA) of a deformable structure depends on the material used (M), the geometry of the deformable zone components (G) and the nature of the deformation force (F D ).Different materials (steel, aluminum alloys, magnesium alloys and composite materials) are used in the body's construction.The impact energy dissipation mechanisms of metallic and composite structures are significantly different.Structural parts of composite materials are brittle and dissipate energy through various combined fracture mechanisms (fiber breakage, matrix cracking, etc.).In contrast, metal structures allow energy dissipation by controlled elastoplastic deformation of the deformation structure components [6][7][8][9].
As mentioned above, the car body consists of a front crumple zone, a passenger safe zone and a rear crumple zone.Due to their required functionality, the material properties of the components used for each zone also differ.For the front and rear crumple zones, dual and multiphase steels (dual-phase steels-DP, complex-phase steels-CP, transformation induced plasticity steels-TRIP, twinning induced plasticity steels-TWIP) are used, which are equipped with very good EA.For the safe zone (safety cage), martensitic steels-MS, hot formed steels-HF, and dual-phase steels >1000 MPa are used, which are characterized by high strength [1,[8][9][10][11].
Dual-phase steels (DP) consist of a fine-grained ferritic matrix in which martensite islands are dispersed.Combining these phases gives the material good fatigue resistance and toughness while maintaining good formability and weldability.DP steels are produced by controlled cooling from the austenitic phase or a two-phase ferritic and austenitic phase.A certain volume fraction of austenite is converted to ferrite, and on sudden cooling, the residual austenite is converted to martensite.The total martensite content in DP steels is in the range of 5-30%.As the proportion of the martensitic phase increases, the strength of the steel increases.High strain hardening exponent n values and high EA of impact energy characterize DP steels [8,12,13].These steels are produced in standardized grades DP350/600, DP 500/800, and DP 600/980.They are used for floor parts, rear bulkheads and deformation stiffeners, front chassis frame, engine mounts, etc. [12].
High-strength multiphase steels with transformation-induced plasticity (TRIP) exhibit an excellent combination of strength and ductility.Their microstructure consists of islands of residual austenite (7-15%) and bainite (30-35%) dispersed in a soft ferritic matrix (50-55%).In some types of TRIP steels, martensite (1-5%) is present.During plastic deformation, austenite is transformed into martensite (TRIP effect).This effect is not concentrated in a particular area but is uniformly distributed throughout the volume.Plates of martensite formed in the slip planes are an obstacle to the movement of dislocations.A martensitic transformation is indicated at its front if a crack develops in the deformed zone during stressing.Plates of martensite then prevent further crack propagation [8,[12][13][14][15][16][17].TRIP steels are produced in the following standardized grades: TRIP 350/600, TRIP 400/700, TRIP 450/800 and TRIP 600/980.Due to their high EA, these steels are used for vehicle body deformable components (transverse and longitudinal beams, B-pillar stiffeners, bumpers, door sills, etc.) [12].
Designers and material engineers frequently face the dilemma of balancing the conflicting demands of vehicle design, weight and crashworthiness.When selecting an appropriate combination of materials for the body deformation structure components, designers need information on their resistance to deformation under compressive or bending forces (Figure 1) and information on their processability (formability, weldability, etc.)-Figure 1 [4,[18][19][20].
This information can be predicted based on the results of the tensile test, three-point bending test, compression test (by squeezing the specimen) and numerical simulation.The present paper aims to draw attention to the possibilities of experimental and numerical methods to predict the stress-strain characteristics of materials such as deformation resistance force-F R or F D , the absorbed energy-EA, the nominal value of deformation resistance force-F MR , the efficiency of deformation resistance force-E FR and the charac-teristic of the stiffness of the deformation structure-CS.The above characteristics were investigated in the research conducted by tensile test, three-point bending test and numerical simulation of three-point bending [4,[21][22][23].This information can be predicted based on the results of the tensile test, three-point bending test, compression test (by squeezing the specimen) and numerical simulation.The present paper aims to draw attention to the possibilities of experimental and numerical methods to predict the stress-strain characteristics of materials such as deformation resistance force-FR or FD, the absorbed energy-EA, the nominal value of deformation resistance force-FMR, the efficiency of deformation resistance force-EFR and the characteristic of the stiffness of the deformation structure-CS.The above characteristics were investigated in the research conducted by tensile test, three-point bending test and numerical simulation of three-point bending.[4,[21][22][23]

Materials
The experimental research of stress-strain characteristics was conducted on sheets of micro-alloyed steel HSLA 240, dual-phase steel DP 600 and multiphase steel TRIP 400/700, which are used in the construction of car bodies.HSLA 240 micro-alloyed steel microstructure is formed by fine equiaxial ferrite grains and small pearlite and tertiary cementite formations-Figure 2a.The microstructure of DP 600 dual-phase steel is a fine-grained polyhedral formed by equiaxial grains of ferrite and very fine grains of martensite-Figure 2b.The TRIP 400/700 microstructure is fine-grained, ferritic, and homogeneous, in which the structural components are arranged in rows.Minority very fine formations of residual austenite and martensite occur in the structure-Figure 2c.
The mechanical properties of the materials were determined by tensile testing on a TIRATEST 2300 testing machine (TIRA Maschinenbau GmbH, Rauenstein, Germany) under the provisions of ISO 6892-1:2019 and ISO 10275:2020.The TIRATES 2300 tensile machine is equipped with a tensometric load cell and extensometer.Five samples were tested in the experiment when the strain rate was 0.002 s −1 .Then, average values of mechanical properties and their standard deviation were calculated from five measurements.Figure 3 shows the dependences of both the engineering stress σE and the true stress σT on the engineering strain εE and the true strain φT, which were obtained from the testing machine recordings.The measured values of the mechanical properties: yield strength Re, tensile

Materials
The experimental research of stress-strain characteristics was conducted on sheets of micro-alloyed steel HSLA 240, dual-phase steel DP 600 and multiphase steel TRIP 400/700, which are used in the construction of car bodies.HSLA 240 micro-alloyed steel microstructure is formed by fine equiaxial ferrite grains and small pearlite and tertiary cementite formations-Figure 2a.The microstructure of DP 600 dual-phase steel is a finegrained polyhedral formed by equiaxial grains of ferrite and very fine grains of martensite-Figure 2b.The TRIP 400/700 microstructure is fine-grained, ferritic, and homogeneous, in which the structural components are arranged in rows.Minority very fine formations of residual austenite and martensite occur in the structure-Figure 2c.
The mechanical properties of the materials were determined by tensile testing on a TIRATEST 2300 testing machine (TIRA Maschinenbau GmbH, Rauenstein, Germany) under the provisions of ISO 6892-1:2019, ISO 10113:2020, and ISO 10275:2020 [24][25][26].The TIRATES 2300 tensile machine is equipped with a tensometric load cell and extensometer.Five samples were tested in the experiment when the strain rate was 0.002 s −1 .Then, average values of mechanical properties and their standard deviation were calculated from five measurements.Figure 3 shows the dependences of both the engineering stress σ E and the true stress σ T on the engineering strain ε E and the true strain φ T , which were obtained from the testing machine recordings.The measured values of the mechanical properties: yield strength R e , tensile strength R m , uniform ductility A U , total ductility A 80 , material constant K 0.1-0.2 and n 0.1-0.2 are presented in      Note: yield strength Re, tensile strength Rm, uniform ductility AU, total ductility A80, Kinterval-strength coefficient, ninterval-strain hardening exponent, 0.002-F-interval within 0.002 strain to fracture strain, 0.05-F-interval within 0.05 strain to fracture strain, 0.1-F-interval within 0.1 strain to fracture strain, 0.1-0.2-interval0.1 within to 0.2 strain, strength coefficient K0.1-0.2, n0.1-0.2 strain hardening exponent, EACTM.t/V0absorbed energy calculated by trapezoidal method according to Equation (2) from the tensile test record from the area under the σT to φT true stress-true strain curve, EACT.t/V0 absorbed energy calculated according to Equation ( 6) from the tensile test record from the area under the σT to φT true stress-true strain curve, V0 unit volume.As the basic deformation characteristic of materials used in body construction, EA expresses the ability to absorb kinetic energy E K into their failure.As it follows from the physio-metallurgical nature of strength and plasticity, EA is closely related to the material's toughness [27,28].Toughness can be characterized as the mechanical work related to the unit volume to failure of the material-EA/V 0 .EA E.t /V 0 can be determined from the tensile test record from the area under the σ E to ε E stress-strain curve-Figure 4, where the dependence of σ E on ε E is presented by the solid line and the dependence of σ T on φ T is presented by the dashed line.Then, EA Eng.t /V 0 was evaluated as follows: Similarly, the EA CTM.t /V 0 characteristics can be calculated by trapezoid method from the dependence of σ T on φ T as follows: As can be observed in Figure 3, the σ T to φ T curve is less flattened than the σ E to ε E dependence curve.This is probably since in the uniform strain area, the cross-section shrinks on one side, but on the other side, there is an intense strain hardening.At the point when σ T starts to decrease, the formation of a neck and, consequently, the failure of the specimen occurs.The strain course in the area of neck formation was imaged with the ARAMIS 3-D optical system.The ARAMIS SRX measuring system with camera resolution of 12 Mpix and 8 GB of the internal memory was used for DIC testing.The maximum sensing speed was 335 Hz a few seconds before specimen fracture, while from the beginning the sensing speed was set to 20 Hz.The sample before the test without the speckles is shown on the bottom, and the sample with speckles on the white background after the fracture is shown on the top-Figure 4.
Similarly, the EACTM.t/V0characteristics can be calculated by trapezoid method from the dependence of σT on φT as follows: As can be observed in Figure 3, the σT to φT curve is less flattened than the σE to εE dependence curve.This is probably since in the uniform strain area, the cross-section shrinks on one side, but on the other side, there is an intense strain hardening.At the point when σT starts to decrease, the formation of a neck and, consequently, the failure of the specimen occurs.The strain course in the area of neck formation was imaged with the ARAMIS 3-D optical system.The ARAMIS SRX measuring system with camera resolution of 12 Mpix and 8 GB of the internal memory was used for DIC testing.The maximum sensing speed was 335 Hz a few seconds before specimen fracture, while from the beginning the sensing speed was set to 20 Hz.The sample before the test without the speckles is shown on the bottom, and the sample with speckles on the white background after the fracture is shown on the top-Figure 4. Figures 5-7 shows a color map of the strain distribution at the point of deformation before failure of the sample.The strong (deep) blue color shows the undeformed area of the sample, the lighter blue color shows the area in which the localization of plastic deformation has not yet occurred, the green color shows the area of the beginning of the localization of plastic deformation, i.e., local narrowing of the cross-section of the sample begins (forming a neck), the yellow color shows the area of the formation of a gradual local narrowing of the cross-section, and the red shows the area just before the failure of the sample.The amount of plastic strain is expressed by the strains φT1, φT2 and φT3 at the limit points SP, WP, CP-Table 2. As can be seen in Figures 5-7, in the blue region, the width of the specimen is the same along the length of the specimen.In the blue region, it can be seen that there is a slight narrowing of the sample width in the deformation localization area, in the yellow region the narrowing of the width is more pronounced than in the green region, and in the red region in the middle of the sample the narrowing is the most pronounced.The strain development at the location of plastic deformation was expressed by the coefficient of normal anisotropy by the Lankford's anisotropy parameter r, which Figures 5-7 shows a color map of the strain distribution at the point of deformation before failure of the sample.The strong (deep) blue color shows the undeformed area of the sample, the lighter blue color shows the area in which the localization of plastic deformation has not yet occurred, the green color shows the area of the beginning of the localization of plastic deformation, i.e., local narrowing of the cross-section of the sample begins (forming a neck), the yellow color shows the area of the formation of a gradual local narrowing of the cross-section, and the red shows the area just before the failure of the sample.The amount of plastic strain is expressed by the strains φ T1 , φ T2 and φ T3 at the limit points SP, WP, CP-Table 2. As can be seen in Figures 5-7, in the blue region, the width of the specimen is the same along the length of the specimen.In the blue region, it can be seen that there is a slight narrowing of the sample width in the deformation localization area, in the yellow region the narrowing of the width is more pronounced than in the green region, and in the red region in the middle of the sample the narrowing is the most pronounced.The strain development at the location of plastic deformation was expressed by the coefficient of normal anisotropy by the Lankford's anisotropy parameter r, which is given by the ratio of the deformation in the direction of the width φ 2 of the sample to the deformation in the direction of the thickness φ 3 of the sample: The r values calculated for individual areas are shown in Table 2.
is given by the ratio of the deformation in the direction of the width φ2 of the sample to the deformation in the direction of the thickness φ3 of the sample: The r values calculated for individual areas are shown in Table 2.If σT from Krupkovsky Equation ( 5) is expressed as follows: Consequently, after substituting σT from Equation ( 5) into Equation ( 4), the following is the result: is given by the ratio of the deformation in the direction of the width φ2 of the sample to the deformation in the direction of the thickness φ3 of the sample: The r values calculated for individual areas are shown in Table 2.The deformation work WD respectively EATC.t can be calculated as the area under the dependence curve σT − φT from the Equation (4): If σT from Krupkovsky Equation ( 5) is expressed as follows: Consequently, after substituting σT from Equation ( 5) into Equation ( 4), the following is the result: is given by the ratio of the deformation in the direction of the width φ2 of the sample to the deformation in the direction of the thickness φ3 of the sample: The r values calculated for individual areas are shown in Table 2.The deformation work WD respectively EATC.t can be calculated as the area under the dependence curve σT − φT from the Equation (4): If σT from Krupkovsky Equation ( 5) is expressed as follows: Consequently, after substituting σT from Equation ( 5) into Equation ( 4), the following is the result: The deformation work W D respectively EA TC.t can be calculated as the area under the dependence curve σ T − φ T from the Equation (4): If σ T from Krupkovsky Equation ( 5) is expressed as follows: Consequently, after substituting σ T from Equation (5) into Equation ( 4), the following is the result: From the measured values of forces and elongations, σ Ti values were calculated using Equation ( 7): Calculating φ Ti strains using the Equation ( 8): where F Ti is a true instantaneous force; φ Ti is the true strain value (φ Ti = ln (1 + ∆L i /L 0 )); ∆L i is the elongation of the specimen under uniaxial tensile stress; L 0 is the specific length of the specimen before deformation; S 0 is the initial specimen cross-section; and A U,i is uniform ductility.The values of K, n obtained by regression analysis at different intervals and the calculated values of EA TRA.t /V 0 using regression analysis under uniaxial tensile stress are given in Table 3. Note: K interaval -strength coefficient; n .interval-strain hardening exponent; EA TRAinterval -energy absorption evaluated by regression analysis; 0.002-F-interval within 0.002 strain to fracture strain; 0.05-F-interval within 0.05 strain to fracture strain; 0.1-F-interval within 0.1 strain to fracture strain.
As mentioned above, in car collisions, the bumper, cross members, sill, B-pillar of the door and roof reinforcement are subjected to bending force.These parts are connected to the body frame structure.In order to make the bending stresses on the deformation components as close as possible to the actual stresses, the ends of the specimens were mechanically clamped in the jaws-Figure 8.For the experimental research, specimens with a width of 30 mm and a length of L = 300 mm were used.Before the deformation of the specimens, a deformation network was etched on their surface.The specimens were deformed to failure at a testing machine traverse speed of 10 mm.min −1 .The bending punch and strain gauge sensor were attached to the moving crossbar of the TIRATEST 2300 testing machine.The machine allows recording the bending force F b and the traverse distance x-Figure 9. Points S0, SZ, WZ and CZ were marked on the curve of the dependence of F b on the strain path x, which delimit the individual deformable zones.F b in the range between points X0 and SZ increases linearly as a function of x pl .The area between these points can be referred to as the safe zone.This zone has a sufficient margin of plastic properties until failure.The curve starts to flatten gradually in the area between the SZ and WZ points, and a neck forms.This area can be described as a warning zone.In the area between the WZ and CP points, the curve of the dependence of F p on x pl is quite flattened.The reserve of plastic properties is exhausted.This area can be referred to as the critical zone.In the area beyond the CP point, the specimen is broken-fracture.The deformation work W D,b or the absorbed energy EA b.CTM was evaluated from the F b record of the x pl test by three-point bending combined stretching the specimen using Equation (9)-Figure 9: whereas F b (x i ) is the average value of the force in the selected interval, x i is the midpoint of the chosen interval [x i−1 ,x i ], and ∆x = x i − x i−1 .
Materials 2024, 17, x FOR PEER REVIEW 9 of 25 on the strain path x, which delimit the individual deformable zones.Fb in the range between points X0 and SZ increases linearly as a function of xpl.The area between these points can be referred to as the safe zone.This zone has a sufficient margin of plastic properties until failure.The curve starts to flatten gradually in the area between the SZ and WZ points, and a neck forms.This area can be described as a warning zone.In the area between the WZ and CP points, the curve of the dependence of Fp on xpl is quite flattened.The reserve of plastic properties is exhausted.This area can be referred to as the critical zone.In the area beyond the CP point, the specimen is broken-fracture.The deformation work WD,b or the absorbed energy EAb.CTM was evaluated from the Fb record of the xpl test by three-point bending combined stretching the specimen using Equation ( 9)-Figure 9: The vehicle body consists of a series of deformable zones [26], each designed to resist deformation until a certain stress level on the part is reached.The deformation elements' resistance force (FR.b = Fb) to plastic deformation is limited by the material properties and geometry.A larger resistance (FR.b) of a deformation element under impact increases the deceleration, which may cause the element to exceed its biomechanical limits.It can be seen from Figure 10 that the area under the curve of the dependence of the bending force Fb on the maximum deformation path xpl.max can be simply calculated as the area of a triangle as follows: A Fbmax can be expressed as a function of the stiffness characteristic: If we substitute Expression (11) for Fbmax in Equation ( 10), we derive the following Equation ( 12): Regression analysis can also determine the stiffness characteristics of CSb,RA.When applying it, the following principles should be observed [27]: • the interval of the strain path xpl in linear regression should be chosen as wide as possible, The vehicle body consists of a series of deformable zones [29], each designed to resist deformation until a certain stress level on the part is reached.The deformation elements' resistance force (F R.b = F b ) to plastic deformation is limited by the material properties and geometry.A larger resistance (F R.b ) of a deformation element under impact increases the deceleration, which may cause the element to exceed its biomechanical limits.It can be seen from Figure 10 that the area under the curve of the dependence of the bending force Fb on the maximum deformation path x pl .max can be simply calculated as the area of a triangle as follows: A F bmax can be expressed as a function of the stiffness characteristic: If we substitute Expression (11) for F bmax in Equation ( 10), we derive the following Equation ( 12): It is possible to determine the stiffness characteristic more precisely if we use EA b.CTM for EA b .F .Then, from Equation ( 12), after adjustment, we arrive at Equation ( 13): Regression analysis can also determine the stiffness characteristics of CS b,RA .When applying it, the following principles should be observed [30]: • the interval of the strain path x pl in linear regression should be chosen as wide as possible,

•
the R 2 value of the variance in regression analysis should be greater than >0.95.
Using the stiffness characteristic CS b,RA determined by regression analysis, Equation ( 14) can be used to calculate EA b,RA /V 0 as follows: The determination of the stiffness characteristics CS b.RA for each material by regression analysis is shown in Figure 9.The measured values of F b = F R.b , x pl , calculated values of stiffness characteristics CS b.RA , CS b.EA and EA b.CTM are presented in Table 4.The failure of the specimens occurred on one side in the center of the unsupported area of the specimen, and on the other side of the unsupported specimen, there was only a narrowing of the specimen cross-section-Figures 10-12.From the deformation grid etched on the specimen, the in-plane φ T1 and φ T2 along the specimen in 11-12 sections were evaluated by the Argus 3D optical system as it is shown in Figures 10-12.The value of the third true strain φT3, i.e., the thickness strain of the specimen, was calculated based on the validity of the volume constancy law according to the Equation ( 15): The Equation ( 15) after adjusting φT3: The measured strain values at the boundaries of each zone are shown in Table 5 Table 5.The strain values measured by the 3D optical system Argus.−0.17 ± 0.01 −0.12 ± 0.01 −0.09 ± 0.01

Numerical Simulation of the Three-Point Bending Test
To predict the stress-strain characteristics of the materials used, a numerical simulation of the bending test with fixed ends was also used in the PAM STAMP 2G V2012 software environment.A 3D model of the actual physical bending tool created in Creo 2.0 was used to set up the numerical simulation-Figure 13a.This was saved in *igs format, and the individual functional parts of the tool were imported into the simulation software.A specimen (strip of sheet metal) with dimensions corresponding to the real specimen (with a length of 120 mm, a width of 30 mm and a thickness of 0.75 mm) was also modelled.After importing them into the PAM STAMP 2G V2012 environment, the different parts of the tool were set to a position that corresponded to the real process.Subsequently, the The value of the third true strain φ T3 , i.e., the thickness strain of the specimen, was calculated based on the validity of the volume constancy law according to the Equation ( 15): The Equation ( 15) after adjusting φ T3 : The measured strain values at the boundaries of each zone are shown in Table 5.

Numerical Simulation of the Three-Point Bending Test
To predict the stress-strain characteristics of the materials used, a numerical simulation of the bending test with fixed ends was also used in the PAM STAMP 2G V2012 software environment.A 3D model of the actual physical bending tool created in Creo 2.0 was used to set up the numerical simulation-Figure 13a.This was saved in *igs format, and the individual functional parts of the tool were imported into the simulation software.A specimen (strip of sheet metal) with dimensions corresponding to the real specimen (with a length of 120 mm, a width of 30 mm and a thickness of 0.75 mm) was also modelled.After importing them into the PAM STAMP 2G V2012 environment, the different parts of the tool were set to a position that corresponded to the real process.Subsequently, the properties such as rigidity, movement, friction, etc. of individual parts of the tool were specified-Figure 13b.The distribution of stresses and strains depends on the behavior of the material.The authors [28,29] describe three models of the elastic-plastic behavior of the material in bending of the cross member: elastic-perfect plastic, elastic-linear plastic, and power-law hardening.The material properties in our numerical simulations were defined by the Hill-48 plasticity condition in combination with the Krupkovsky isotropic hardening model (Equation ( 5)).The values of Young's modulus of elasticity E, Poisson's number µ, strength constant K, strain hardening exponent n and the plastic strain ratios r were specified in the material card, and the boundary conditions (friction coefficient, loading direction, holding force, punch speed, etc.) were set in the process tree of the simulation software [30].Surface blank was meshed with square elements of 0.5 mm size and the mesh was checked with no errors.
Measured values of materials-strength constant K, strain hardening exponent n and Lankford's anisotropy parameter r (Table 4)-were set as input values in the numerical simulation of three-point bending.Subsequently, the boundary conditions corresponding to the conditions of the three-point bending test were set.Before the true strain, the plate strip (specimen) was held at the edge by a force (300 kN) to prevent it from pulling out of the jaw.Subsequently, a bending path of 120 mm was defined.The deformation of the specimen was followed by the actual calculation of the stress-strain characteristics (Fb,xpl, in-plane strains φ1T, φ2T and the thickness strain φ3T) in 11 steps.Figure 14 to Figure 16 show the dependences of the bending force Fs.b on the punch path xpl.s, which were obtained by numerical simulation.The elastic xel and plastic xpl phases of the strain path are indicated in these figures.When evaluating EAS.b, the stress-strain characteristics were shifted to the left by xel (approximately 8 to 10 mm) i.e., by the restitution work (elastic work EAel, etc.)-Figures 14-16.The distribution of stresses and strains depends on the behavior of the material.The authors [31,32] describe three models of the elastic-plastic behavior of the material in bending of the cross member: elastic-perfect plastic, elastic-linear plastic, and power-law hardening.The material properties in our numerical simulations were defined by the Hill-48 plasticity condition in combination with the Krupkovsky isotropic hardening model (Equation ( 5)).The values of Young's modulus of elasticity E, Poisson's number µ, strength constant K, strain hardening exponent n and the plastic strain ratios r were specified in the material card, and the boundary conditions (friction coefficient, loading direction, holding force, punch speed, etc.) were set in the process tree of the simulation software [33].Surface blank was meshed with square elements of 0.5 mm size and the mesh was checked with no errors.
Measured values of materials-strength constant K, strain hardening exponent n and Lankford's anisotropy parameter r (Table 4)-were set as input values in the numerical simulation of three-point bending.Subsequently, the boundary conditions corresponding to the conditions of the three-point bending test were set.Before the true strain, the plate strip (specimen) was held at the edge by a force (300 kN) to prevent it from pulling out of the jaw.Subsequently, a bending path of 120 mm was defined.The deformation of the specimen was followed by the actual calculation of the stress-strain characteristics (F b,xpl , in-plane strains φ 1T , φ 2T and the thickness strain φ 3T ) in 11 steps.Figures 14-16 show the dependences of the bending force F s.b on the punch path x pl.s , which were obtained by numerical simulation.The elastic x el and plastic x pl phases of the strain path are indicated in these figures.When evaluating EA S.b , the stress-strain characteristics were shifted to the left by x el (approximately 8 to 10 mm) i.e., by the restitution work (elastic work EA el , etc.)-Figures 14-16.Similar to the above, the deformable zones were marked on the curve of the dependence of the deformation force Fb on the strain path xpl.The dependence of Fs.b on xpl.s is linear in the area between the CP and SP points.This zone can be referred to as the safe zone.In this zone, it is possible to fully exploit the properties of metallic materials in controlling the deformation process of body parts in impact.In the area between the SP and WP points, the dependence of Fs.b on xpl.s begins to flatten and the plastic deformation gradually localizes.The specimen at the point of localization begins to narrow.This zone can be referred to as the warning zone (i.e., edge zone), in which the plastic deformation process of the deformed body parts can also still be controlled.In the critical zone between the WP and CP points, a more pronounced flattening of the Fs.b vs, xpl.s dependence occurs.At the CP point, the plastic properties are exhausted; hence, the possibility of controlling the deformation process of body parts during impact is lost.Table 6 shows the measured values of xpl.s, Fs.b, EAs.b, CSs.b in the plastic deformation area.In evaluating these characteristics, the same procedure was used to evaluate the experimentally determined stressstrain characteristics of the materials used.Similar to the above, the deformable zones were marked on the curve of the dependence of the deformation force F b on the strain path x pl .The dependence of F s.b on x pl.s is linear in the area between the CP and SP points.This zone can be referred to as the safe zone.In this zone, it is possible to fully exploit the properties of metallic materials in controlling the deformation process of body parts in impact.In the area between the SP and WP points, the dependence of F s.b on x pl.s begins to flatten and the plastic deformation gradually localizes.The specimen at the point of localization begins to narrow.This zone can be referred to as the warning zone (i.e., edge zone), in which the plastic deformation process of the deformed body parts can also still be controlled.In the critical zone between the WP and CP points, a more pronounced flattening of the F s.b vs, x pl.s dependence occurs.At the CP point, the plastic properties are exhausted; hence, the possibility of controlling the deformation process of body parts during impact is lost.Table 6 shows the measured values of x pl.s , F s.b , EA s.b , CS s.b in the plastic deformation area.In evaluating these characteristics, the same procedure was used to evaluate the experimentally determined stress-strain characteristics of the materials used.

Results and Discussion
The basic deformation characteristic of the materials used in the construction of the body is the ability to absorb (EA) kinetic energy (EK) upon impact to material failure.As follows from Equation (1), the lower predictive ability of EA Eng.t is related to the fact that the initial S 0 cross-section and not the actual S T cross-section is used in its determination.The predictive ability of the dependence σ T on φ T using the Krupkovsky equation was also verified outside the interval recommended by ISO 10275:2020 [26], namely in the following strain intervals: φ T from 0.002 to the uniform elongation A U (pre-break φ T,F ), φ T,F from 0.5 to the uniform elongation A U (pre-break φ T,F ), and φ T from 0.1 to the uniform elongation A U (fracture or pre-break φ T,F )-Figures 17-19.

Results and Discussion
The basic deformation characteristic of the materials used in the construction of the body is the ability to absorb (EA) kinetic energy (EK) upon impact to material failure.As follows from Equation (1), the lower predictive ability of EAEng.t is related to the fact that the initial S0 cross-section and not the actual ST cross-section is used in its determination.The predictive ability of the dependence σT on φT using the Krupkovsky equation was also verified outside the interval recommended by ISO 10275:2020, namely in the following strain intervals: φT from 0.002 to the uniform elongation AU (pre-break φT,F), φT,F from 0.5 to the uniform elongation AU (pre-break φT,F), and φT from 0.1 to the uniform elongation AU (fracture or pre-break φT,F)-Figures 17-19.

Results and Discussion
The basic deformation characteristic of the materials used in the construction of the body is the ability to absorb (EA) kinetic energy (EK) upon impact to material failure.As follows from Equation (1), the lower predictive ability of EAEng.t is related to the fact that the initial S0 cross-section and not the actual ST cross-section is used in its determination.The predictive ability of the dependence σT on φT using the Krupkovsky equation was also verified outside the interval recommended by ISO 10275:2020, namely in the following strain intervals: φT from 0.002 to the uniform elongation AU (pre-break φT,F), φT,F from 0.5 to the uniform elongation AU (pre-break φT,F), and φT from 0.1 to the uniform elongation AU (fracture or pre-break φT,F)-Figures 17-19.The designer's goal during the process of car body design is to reduce its weight and increase resistance to impact.Criteria for optimization for the impact resistance of body components are energy absorption, stiffness, etc.The absorption capacity EACTM.t and EAb.CTM of the investigated materials was determined from the record of the tensile test and the three-point bending test using the trapezoidal method using Equations ( 2) and ( 9).The results obtained by the trapezoidal method were considered as a reference when comparing several assessment methods.The Figure 20 shows a comparison of the absorption capacity of micro-alloyed steel HSLA 240, two-phase steel DP 600, and multiphase steel TRIP 400/700 from the tensile test record.From the comparison of EACTM.tDP 600 against micro-alloyed steel HSLA 240, it can be seen that this characteristic has increased by 30% and compared to TRIP 400/700 by up to 79%.The increase in absorption capacity of DP 600 and TRIP 400/700 steel compared to The designer's goal during the process of car body design is to reduce its weight and increase resistance to impact.Criteria for optimization for the impact resistance of body components are energy absorption, stiffness, etc.The absorption capacity EA CTM.t and EA b.CTM of the investigated materials was determined from the record of the tensile test and the three-point bending test using the trapezoidal method using Equations ( 2) and ( 9).The results obtained by the trapezoidal method were considered as a reference when comparing several assessment methods.The Figure 20 shows a comparison of the absorption capacity of micro-alloyed steel HSLA 240, two-phase steel DP 600, and multiphase steel TRIP 400/700 from the tensile test record.The designer's goal during the process of car body design is to reduce its weight and increase resistance to impact.Criteria for optimization for the impact resistance of body components are energy absorption, stiffness, etc.The absorption capacity EACTM.t and EAb.CTM of the investigated materials was determined from the record of the tensile test and the three-point bending test using the trapezoidal method using Equations ( 2) and ( 9).The results obtained by the trapezoidal method were considered as a reference when comparing several assessment methods.The Figure 20 shows a comparison of the absorption capacity of micro-alloyed steel HSLA 240, two-phase steel DP 600, and multiphase steel TRIP 400/700 from the tensile test record.From the comparison of EACTM.tDP 600 against micro-alloyed steel HSLA 240, it can be seen that this characteristic has increased by 30% and compared to TRIP 400/700 by up to 79%.The increase in absorption capacity of DP 600 and TRIP 400/700 steel compared to From the comparison of EA CTM.t DP 600 against micro-alloyed steel HSLA 240, it can be seen that this characteristic has increased by 30% and compared to TRIP 400/700 by up to 79%.The increase in absorption capacity of DP 600 and TRIP 400/700 steel compared to micro-alloyed steel is related to the strength coefficient K and strain hardening exponent n.In the case of DP 600, the value of K was greater by 47% and in the case of n by 13%.In the case of TRIP 400/700, the value of K was greater by up to 204% and n was greater by 61%.From the mentioned comparison, it follows that materials with a two-and multiphase structure strengthen more intensively than materials with a single-phase ferritic structure.The evaluation of EA CTM.t by the trapezoidal method is accurate but time-consuming, therefore the Relation (6) was derived.From the comparison of EA CTM.t and EA CT.t values calculated according to Relation (6) shows that agreement between these characteristics was achieved in the range from 1% to 6%.When qualitatively evaluating the dependence model σ T on φ T on Figure 21, we observe that R 2 for the dependence of EA RA.t /V 0 and the measured values of EA CTM.t /V 0 takes the value of 0.9998.Similarly with the dependence of EA CT.t /V 0 on EA CTM.t R 2 was 0.9932.The dependence of EA CT.t /V 0 on EA CTM.t shows a tendency that can be expressed as follows: If this tendency is confirmed or, if it could be, refined, then after substituting for EA CT.t /V 0 the expression from the right side of Equation ( 6) we derive: Deformable parts of the body (bumper, B-pillar, door and roof reinforcements, etc.) are exposed to a bending force during an impact.In this case, the results of EA b.CTM evaluated by the trapezoidal method were considered as a reference when comparing.The absorption capacity (EA b.CTM ) of the materials used was evaluated from the recording of the dependence of the bending force (F b ) on the deformation path x.The Figure 22 shows a comparison of EA b.CTM of the investigated materials.From the comparison, it can be seen that the material DP 600 shows 1.3 times greater absorption potential and the material TRIP 400/700 approx.1.9 times greater impact energy absorption potential than HSLA 240 material.
Materials 2024, 17, x FOR PEER REVIEW 18 micro-alloyed steel is related to the strength coefficient K and strain hardening exp n.In the case of DP 600, the value of K was greater by 47% and in the case of n by 13 the case of TRIP 400/700, the value of K was greater by up to 204% and n was great 61%.From the mentioned comparison, it follows that materials with a two-and m phase structure strengthen more intensively than materials with a single-phase fe structure.The evaluation of EACTM.t by the trapezoidal method is accurate but time suming, therefore the Relation (6) was derived.From the comparison of EACTM.t and values calculated according to Relation (6) shows that agreement between these ch teristics was achieved in the range from 1% to 6%.When qualitatively evaluating th pendence model σT on φT on Figure 21 Deformable parts of the body (bumper, B-pillar, door and roof reinforcements are exposed to a bending force during an impact.In this case, the results of EAb.CTM e ated by the trapezoidal method were considered as a reference when comparing.Th sorption capacity (EAb.CTM) of the materials used was evaluated from the recording o dependence of the bending force (Fb) on the deformation path x.The Figure 22 sho  Another deformation characteristic is stiffness (CS).The CSb.EA value calculated according to Relation (13) was considered as a reference value when comparing the stiffness of the investigated materials.It can be seen from the CSb.EA comparison shown in the Figure 23, that DP 600 material exhibits 1.54 times greater stiffness and TRIP 400/700 approximately 1.74 times greater stiffness than HSLA 240.Based on the above, it can be assumed that DP 600 materials and TRIP 400/700 have the potential to reduce the weight of body parts (using thinner wall thickness components) while maintaining or even improving impact resistance properties compared to HSLA 240.Furthermore, the possibility of using regression analysis in determining the characteristics of CSb.RA was verified.A comparison of CSb.RA and CSb.EA characteristics show that the agreement between these characteristics ranges from 86% to 95%.The possibility of Another deformation characteristic is stiffness (CS).The CS b.EA value calculated according to Relation (13) was considered as a reference value when comparing the stiffness of the investigated materials.It can be seen from the CS b.EA comparison shown in the Figure 23, that DP 600 material exhibits 1.54 times greater stiffness and TRIP 400/700 approximately 1.74 times greater stiffness than HSLA 240.Based on the above, it can be assumed that DP 600 materials and TRIP 400/700 have the potential to reduce the weight of body parts (using thinner wall thickness components) while maintaining or even improving impact resistance properties compared to HSLA 240.Another deformation characteristic is stiffness (CS).The CSb.EA value calculated according to Relation (13) was considered as a reference value when comparing the stiffness of the investigated materials.It can be seen from the CSb.EA comparison shown in the Figure 23, that DP 600 material exhibits 1.54 times greater stiffness and TRIP 400/700 approximately 1.74 times greater stiffness than HSLA 240.Based on the above, it can be assumed that DP 600 materials and TRIP 400/700 have the potential to reduce the weight of body parts (using thinner wall thickness components) while maintaining or even improving impact resistance properties compared to HSLA 240.Furthermore, the possibility of using regression analysis in determining the characteristics of CSb.RA was verified.A comparison of CSb.RA and CSb.EA characteristics show that the agreement between these characteristics ranges from 86% to 95%.The possibility of Furthermore, the possibility of using regression analysis in determining the characteristics of CS b.RA was verified.A comparison of CS b.RA and CS b.EA characteristics show that the agreement between these characteristics ranges from 86% to 95%.The possibility of predicting EA and CS was also verified based on the results of numerical simulation.For EA, concordance with EA s.b ranged from 86% to 90%, and for CS .s.b ranged from 64% to 67%.The issue of verifying the results of numerical simulation will be addressed in further research.
The characteristics of EA and CS can guide designers when selecting materials for deformable body parts.For example, if F b loads the deformation element, the deceleration can be expressed from the energy balance condition: Assuming that a vehicle of weight m moves on a path x pl with acceleration a acquires kinetic energy: This E K on impact is consumed by the W D or absorbed by the deforming structure in the x pl path.The vehicle moves on this path with a deceleration a.When substituted into the energy equilibrium condition (Equation ( 19), we obtain: After adjusting the deceleration a is as follows: The effect of material substitution on deceleration can be expressed by the deceleration change index Ia: whereas a i is a deceleration in the use of innovative material, a ref is the deceleration when using the reference material.
If we assume that on the same vehicle (of the same weight), the reference material of micro-alloyed HSLA 240 steel will be replaced with DP 600 dual-phase steel, then after substituting for a i in Equation ( 23) and adjusting, the result is as follows: The comparison of the I a indices calculated using the CS b.EA Equation ( 13) and the regression analysis CS b.RA indicates that the differences are very small, ranging from −1% to 3%.When substituting HSLA 240 with DP 600, a change in the deceleration index of up to 30% can be expected, and when substituting HSLA 240 with TRIP 400/700, a change in deceleration of up to 69% can be expected.
Another important characteristic is the strain path x pl , at which controlled deformation of the deformation structure occurs.From the comparison of the experimentally measured x pl values and the calculated x pl,s shown in Table 6, it can be observed that there was a difference of 4% for the HSLA 240 material with a single-phase structure and a difference of 27% to 32% for the DP 600 material with a two-phase structure.The largest difference between x pl and x pl,s was observed for the TRIP 400/700 material with a three-phase structure (30% to 43%)-Table 6.We assume that these differences are related to the model used for prediction in the program file.Although plastic deformation is localized in the warning zone between the SP and WP points (the beginning of the neck formation was observed on the specimen), it is still possible to control the body parts' deformation processes at impact effectively.
For predicting deformation states can be used FLC [34][35][36].These curves allow the separation of the different transition zones from the safe deformation behavior of the material (with a certain margin) to the states where plastic deformation is localized, i.e., neck formation and subsequent failure.In FLC design, it is important to know the value of FLC0 i.e., the state when at φ 2T = 0 the value of φ 1T (FLC0 = φ 1T ) is determined.The example is shown in Figure 24 for DP 600 material.Based on the results of [34], the value of FLC0 for HSLA 260 material is about 0.3, based on [35], the value of FLC0 for DP 600 material is about 0.2 and based on [36] for TRIP 400/600 material the value is about 0.2.These FLC0 CP values for separating the critical zone from the warning zone (FLC0 WP ) are recommended to be shifted by 10%-Table 7.  [34][35][36] 0.27 [34] 0.18 [35] 0.18 [36] Materials 2024, 17, x FOR PEER REVIEW 21 of 25 Another important characteristic is the strain path xpl, at which controlled deformation of the deformation structure occurs.From the comparison of the experimentally measured xpl values and the calculated xpl,s shown in Table 6, it can be observed that there was a difference of 4% for the HSLA 240 material with a single-phase structure and a difference of 27% to 32% for the DP 600 material with a two-phase structure.The largest difference between xpl and xpl,s was observed for the TRIP 400/700 material with a threephase structure (30% to 43%)-Table 6.We assume that these differences are related to the model used for prediction in the program file.Although plastic deformation is localized in the warning zone between the SP and WP points (the beginning of the neck formation was observed on the specimen), it is still possible to control the body parts' deformation processes at impact effectively.
For predicting deformation states can be used FLC [31][32][33].These curves allow the separation of the different transition zones from the safe deformation behavior of the material (with a certain margin) to the states where plastic deformation is localized, i.e., neck formation and subsequent failure.In FLC design, it is important to know the value of FLC0 i.e., the state when at φ2T = 0 the value of φ1T (FLC0 = φ1T) is determined.The example is shown in Figure 24 for DP 600 material.Based on the results of [31], the value of FLC0 for HSLA 260 material is about 0.3, based on [32], the value of FLC0 for DP 600 material is about 0.2 and based on [33] for TRIP 400/600 material the value is about 0.2.These FLC0CP values for separating the critical zone from the warning zone (FLC0WP) are recommended to be shifted by 10%-Table 7.  [31][32][33] 0.27 [31] 0.18 [32] 0.18 [33] The curve that separates the critical zone from the warning zone is indicated in the FLD by a dashed line-Figure 24.A comparison of the numerical simulation results presented in Table 7 with the experimental results published in [29][30][31] shows for the HSLA 240 material, with the assumed strain φ2WP = 0, the FLD corresponds to FLC0WP.s= 0. 27; for the DP 600 material the value FLC0WP.s= 0.18 [30]; and for the TRIP 400/700 material the value FLC0WP.s= 0.18 [31].The above comparison shows that the differences of FLC0WP.sobtained by simulation for the DP 600 and TRIP 400/700 material are 50% as large as the data published in [32,33].The curve that separates the critical zone from the warning zone is indicated in the FLD by a dashed line-Figure 24.A comparison of the numerical simulation results presented in Table 7 with the experimental results published in [32][33][34] shows for the HSLA 240 material, with the assumed strain φ 2WP = 0, the FLD corresponds to FLC0 WP.s = 0. 27; for the DP 600 material the value FLC0 WP.s = 0.18 [33]; and for the TRIP 400/700 material the

Materials 2024 , 25 Figure 1 .
Figure 1.Scheme of the three-point bending test and compression test.

Figure 1 .
Figure 1.Scheme of the three-point bending test and compression test.

Figure 3 .
Figure 3. Dependence of engineering σ E and true σ T stress on true φ T strain.

Figure 4 .
Figure 4.The sample before and after the tensile test when recorded by ARAMIS SRX optical system.

Figure 4 .
Figure 4.The sample before and after the tensile test when recorded by ARAMIS SRX optical system.

Figure 5 .
Figure 5. Deformation record by the ARAMIS system along the HSLA 240 specimen during the tensile test.

Figure 6 .
Figure 6.Deformation record by the ARAMIS system along the DP 600 specimen during the tensile test.

Figure 7 .
Figure 7. Deformation record by the ARAMIS system along the TRIP 400/700 specimen during the tensile test.The deformation work WD respectively EATC.t can be calculated as the area under the dependence curve σT − φT from the Equation (4):

Figure 5 .
Figure 5. Deformation record by the ARAMIS system along the HSLA 240 specimen during the tensile test.

Figure 5 .
Figure 5. Deformation record by the ARAMIS system along the HSLA 240 specimen during the tensile test.

Figure 6 .
Figure 6.Deformation record by the ARAMIS system along the DP 600 specimen during the tensile test.

Figure 7 .
Figure 7. Deformation record by the ARAMIS system along the TRIP 400/700 specimen during the tensile test.

Figure 6 .
Figure 6.Deformation record by the ARAMIS system along the DP 600 specimen during the tensile test.

Figure 5 .
Figure 5. Deformation record by the ARAMIS system along the HSLA 240 specimen during the tensile test.

Figure 6 .
Figure 6.Deformation record by the ARAMIS system along the DP 600 specimen during the tensile test.

Figure 7 .
Figure 7. Deformation record by the ARAMIS system along the TRIP 400/700 specimen during the tensile test.

Figure 7 .
Figure 7. Deformation record by the ARAMIS system along the TRIP 400/700 specimen during the tensile test.

Figure 8 .
Figure 8. Principle scheme of the three-point bending test.

Figure 9 .
Figure 9. Record of the bending force on the punch path.

Figure 10 .
Figure 10.Record from the strain evaluation by the 3D optical system Argus-material HSLA 240.Note: The measured strain φ1 values are indicated by a thin solid black line and φ3 value are indicated by a thin solid blue line.

Figure 11 .
Figure 11.Record from the strain evaluation by the 3D optical system Argus-material DP 600.Note: The measured strain φ1 values are indicated by a thin solid black line and φ3 value are indicated by a thin solid blue line.

Figure 10 .
Figure 10.Record from the strain evaluation by the 3D optical system Argus-material HSLA 240.Note: The measured strain φ 1 values are indicated by a thin solid black line and φ 3 value are indicated by a thin solid blue line.

Figure 10 .
Figure 10.Record from the strain evaluation by the 3D optical system Argus-material HSLA 240.Note: The measured strain φ1 values are indicated by a thin solid black line and φ3 value are indicated by a thin solid blue line.

Figure 11 .
Figure 11.Record from the strain evaluation by the 3D optical system Argus-material DP 600.Note: The measured strain φ1 values are indicated by a thin solid black line and φ3 value are indicated by a thin solid blue line.

Figure 11 .
Figure 11.Record from the strain evaluation by the 3D optical system Argus-material DP 600.Note: The measured strain φ 1 values are indicated by a thin solid black line and φ 3 value are indicated by a thin solid blue line.

Figure 12 .
Figure 12.Record from the strain evaluation by the 3D optical system Argus-material TRIP 400/700.Note: The measured strain φ1 values are indicated by a thin solid black line and φ3 value are indicated by a thin solid blue line.

Figure 12 .
Figure 12.Record from the strain evaluation by the 3D optical system Argus-material TRIP 400/700.Note: The measured strain φ 1 values are indicated by a thin solid black line and φ 3 value are indicated by a thin solid blue line.

Materials 2024 ,Figure 13 .
Figure 13.Three-point bending test with fixed ends: (a) 3-D model created in Creo 2.0; (b) simulation model created in Pam Stamp software.

Figure 13 .
Figure 13.Three-point bending test with fixed ends: (a) 3-D model created in Creo 2.0; (b) simulation model created in Pam Stamp software.

Figure 14 .
Figure 14.Comparison of simulation and experiment results for material HSLA 240.

Figure 15 .
Figure 15.Comparison of simulation and experiment results for material DP 600.

Figure 14 . 25 Figure 14 .
Figure 14.Comparison of simulation and experiment results for material HSLA 240.

Figure 15 .
Figure 15.Comparison of simulation and experiment results for material DP 600.

Figure 15 .
Figure 15.Comparison of simulation and experiment results for material DP 600.

Figure 16 .
Figure 16.Comparison of simulation and experiment results for material TRIP 400/700.

Figure 16 .
Figure 16.Comparison of simulation and experiment results for material TRIP 400/700.

Figure 17 .
Figure 17.Dependence of the true stress on the true strain within the strain interval <0.002; φT,F>.

Figure 18 .
Figure 18.Dependence of the true stress on the true strain within the strain interval <0.05; φT,F>.

Figure 17 .
Figure 17.Dependence of the true stress on the true strain within the strain interval <0.002; φ T,F >.

Figure 17 .
Figure 17.Dependence of the true stress on the true strain within the strain interval <0.002; φT,F>.

Figure 18 .
Figure 18.Dependence of the true stress on the true strain within the strain interval <0.05; φT,F>.Figure 18. Dependence of the true stress on the true strain within the strain interval <0.05; φ T,F >.

Figure 18 .
Figure 18.Dependence of the true stress on the true strain within the strain interval <0.05; φT,F>.Figure 18. Dependence of the true stress on the true strain within the strain interval <0.05; φ T,F >.

Figure 19 .
Figure 19.Dependence of the true stress on the true strain within the strain interval <0.1; φT,F>.

Figure 20 .
Figure 20.Comparison of the absorption energy of the investigated materials on the basis of applied models-tensile test.

Figure 19 .
Figure 19.Dependence of the true stress on the true strain within the strain interval <0.1; φ T,F >.

Materials 2024 , 25 Figure 19 .
Figure 19.Dependence of the true stress on the true strain within the strain interval <0.1; φT,F>.

Figure 20 .
Figure 20.Comparison of the absorption energy of the investigated materials on the basis of applied models-tensile test.

Figure 20 .
Figure 20.Comparison of the absorption energy of the investigated materials on the basis of applied models-tensile test.

Figure 21 .
Figure 21.Dependence of the energy absorbing calculated from the regression model, accord Equation (6) models and according to Equation (1) on the energy absorbing according to Equ (2).

Figure 21 .
Figure 21.Dependence of the energy absorbing calculated from the regression model, according to Equation (6) models and according to Equation (1) on the energy absorbing according to Equation (2).

Figure 22 .
Figure 22.Comparison of the absorption energy of the investigated materials on the basis of applied models-three-point bending test.

Figure 23 .
Figure 23.Comparison of the stiffness characteristics of the investigated materials based on the applied models-three-point bending test.

Figure 22 .
Figure 22.Comparison of the absorption energy of the investigated materials on the basis of applied models-three-point bending test.

Materials 2024 , 25 Figure 22 .
Figure 22.Comparison of the absorption energy of the investigated materials on the basis of applied models-three-point bending test.

Figure 23 .
Figure 23.Comparison of the stiffness characteristics of the investigated materials based on the applied models-three-point bending test.

Figure 23 .
Figure 23.Comparison of the stiffness characteristics of the investigated materials based on the applied models-three-point bending test.

Table 1 .
Measured values of mechanical properties of the materials used.

Table 1 .
Measured values of mechanical properties of the materials used.

Table 2 .
The strain values measured with the Aramis optical system.

Table 3 .
Measured values of n, K and energy absorbing at different intervals of true strain.

Table 4 .
Measured and calculated values of the deformation characteristics by the three-point bending test.

Table 6 .
Deformation characteristics determined by simulation of the three-point bending test.

Table 6 .
Deformation characteristics determined by simulation of the three-point bending test.